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Abstract. 

We calculate the k = 2 string tension in SU(4) and SU(5) gauge theories in 3+1 
dimensions, and compare it to the k — 1 fundamental string tension. We find, from 
the continuum extrapolation of our lattice calculations, that 0*^=2/ °~f — 1.40 ± 0.08 in 
the SU(4) gauge theory, and that o~k=2i 1 °f = 1.56±0.10 in SU(5). We remark upon the 
way this might constrain the dynamics of confinement and the intriguing implications 
it might have for the mass spectrum of SU(iV) gauge theories. We also note that these 
results agree closely with the MQ CD-inspired conjecture that the SU(iV) string tension 
varies as oc sin(7rfc/iV). 



1 Introduction 



It is well known that confinement in four dimensional SU(iV) gauge theories leads to 
a linear potential between static test charges in the fundamental representation, so 
that there is a non-zero fundamental string tension af ^ 0. For SU(iV > 4) there is, 
however, the possibility of new stable strings which join test charges in representations 
higher than the fundamental. One can conveniently label these by the way the test 
charge transforms under a gauge transformation, z, that belongs to the centre, Z^, of 
the group. If the charge acquires a factor z k we shall refer to the string as having N- 
ality k, and the corresponding string tension will be denoted by cr^. So the fundamental 
string has iV-ality k = 1 and we shall refer to it as 07 or Ofc=ij interchangeably. 

Because gluons transform trivially under the centre, such a fc-string will not be 
screened down to a ^'-string, if k' 7^ k. If, however, z k = z k Vz then the A;-string can 
transform into a fc'-string and so one trivially has that a^ = ay. For example, in the 
case of SU(3) a k = 2 string joining a diquark source to a distant anti-diquark source 
will be simply the usual fundamental k — 1 string over almost all of the separation. 
(In SU(3) three fundamental strings can form a colour singlet.) We have to go to at 
least SU(4) in order to have the possibility of a genuinely different k = 2 string, and 
to SU(6) for a k = 3 string. 

One can form a fc-string by simply using static test charges that consist of k fun- 
damental charges. So the most trivial possibility is that a '^-string' is nothing more 
than k separate fundamental strings joining the two test charges. In this case we will 
have Ofe = kerf. So for example, if <Jk=2 ^2<r/ then the k = 2 string is just like two 
k = 1 strings and we have nothing new. There is however a non-trivial possibility. If 
<7fc = 2 < 2a f then the k = 2 string will be stable and will constitute a new kind of string. 
We may regard it as a bound state of two k = 1 strings. Similar comments apply to 
higher values of k. 

As emphasised in [|lj, one reason why the lattice calculation of a^ in SU(iV) gauge 
theories is particularly interesting is that it tests the conjecture || that the ratio a^jaf 
has a universal value in 'QCD-like' theories: 

Ok = sin% 

a f sinf KJ 

This conjecture is based on calculations in brane (M-)theory Jl||. The resulting variants 
of QCD go under the generic name of MQCD. One finds that in MQCD with M = 1 
supersymmetry || |3|, the string tensions of fc-strings satisfy eqn(|lj). The same is true 
for M = 2, for Af = 2 softly broken to Af = 1 and for non-supersymmetric MQCD 
[|]. This is so even though the actual values of at differ in these theories. Hence the 
conjecture that the ratio of string tensions is 'universal' in QCD-like theories and 
takes the value in eqn ([!]). 

Another reason for being interested in the value of a^ is that the existence of k- 
strings may have a striking impact on the mass spectrum of SU(iV) gauge theories 
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in the large- N limit [[|. A natural model for the glueball states that make up the 
spectrum of pure gauge theories is as quantised closed strings of fundamental flux. 
Such a picture is most compelling for highly excited states which are large compared 
to the intrinsic width of the flux tube, but may be a reasonably accurate representation 
for the low-lying spectrum as well. A simple and explicit model based on this idea was 
suggested some time ago || and comparisons || |7|, |8| with the detailed SU(iV) lattice 
mass spectra available in 2+1 dimensions show that it works reasonably well. As 
emphasised in || one can build towers of glueball states on any stable fc-string and 
these will be identical up to a rescaling of the overall mass scale by a factor of Jcr^/af 
(if one neglects the finite width of the flux tube). This argument neglects mixing and 
decay, but this becomes appropriate as N — > oo. Such a mass spectrum would provide 
remarkable evidence for the underlying string structure of the glueball states. All this 
does however depend on the existence of stable fc-strings and the details of the SU(iV) 
spectrum will depend on the actual values of the corresponding string tensions. 

Finally we remark that knowing the value of Ofc/07 will constrain the details of any 
proposed confinement mechanism. 

In this paper we will present some results on the k = 2 string tension and will 
compare it to the fundamental (k = 1) tension. We do this for both SU(4) and SU(5), 
although our calculations will be more precise in the former case. There exist previous 
calculations |I| of the lattice SU(4) k = 2 string tension which show that < 2a y. 
However these calculations were performed at quite a high temperature (although still 
in the confining phase) and were not accurate enough to make a useful statement about 
the continuum limit. We intend to perform calculations that are accurate enough, in 
the continuum limit, to provide a significant test of the MQCD conjecture in eqn(|l]). 
These calculations are part of our ongoing study [|Il]] of the physical properties of 



D=3+l SU(./V) gauge theories, and a similar calculation is in progress elsewhere [12 



2 Calculations 

To calculate we first construct a fc-string that closes upon itself through a spatial 
boundary (so no explicit static test charges are needed). We then calculate the corre- 
lation function of two such fc-strings as a function of their (Euclidean) time separation, 
using standard Monte Carlo simulation methods. At sufficiently large separation t, the 
correlation function will fall as a simple exponential oc exp(— m^t) where mk is the mass 
of the lightest /c-string state. Linear confinement means that this mass will be linear 
in the length of the string, which we denote by / = aL where a is the lattice spacing 
and L is the lattice spatial size in lattice units. In addition to this linear dependence 
the first correction is universal ]n| and we assume that the universality class is that 
of a simple bosonic string. (There is evidence from many previous SU(2) and SU(3) 
lattice calculations that points to this.) Using the string correction appropriate to our 
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closed periodic string [14| we therefore have 



n 7T 

am k = a a k L - — + . . . (2) 

where everything has been expressed in lattice units. So once we have calculated the 
mass amk we can use eqn(^) to extract a value for a 2 0)~. 

In practice a calculation using the simplest lattice fc-string operator is inefficient 
because the overlap onto the lightest string state is small and so one has to go to large 
values of t before the contribution of excited states has died away; and there the signal 
disappears into the statistical noise. There are standard methods ]K| for curing this 



problem, using blocked (smeared) link operators and variational techniques. This is 
described in detail in, for example, || for the case of a k — 1 string. Let Trl\ be such 
a string, so that l\ is the string before one closes it by taking the trace. Then we can 
form trial k = 2 strings from (Tr^) 2 and Tr (ZiZi) . In the variational calculation we 
allow arbitrary linear combinations of these in the search for the lightest state. (In 
practice we expect this to belong to the appropriate representation of SU(iV).) Details 
of this, and of our choice of Monte Carlo algorithm, will be given elsewhere [|TT] . 



Our SU(4) calculations are summarised in Table [T| They span a range of couplings 
that corresponds roughly, to the range (3 G [5.7,6.0] in the case of SU(3) |T7|]. This 
is well beyond the bulk phase transition due to the non-trivial phase structure in the 
fundamental-adjoint coupling plane This phase structure will be more complex 
as iV increases, since more couplings, appropriate to representations other than just 
the fundamental and adjoint, become relevant, but its effects are expected to remain 
confined to the region between weak and strong coupling. (The phase structure is 
believed to be due to the differing sensitivity of plaquettes in different representations 
to the presence of Z%. monopoles of various k, and there are more possibilities as N 
increases.) 

The volumes used in these calculations have been chosen so that if we were working 
in SU(3) on the same volumes (when expressed in physical units e.g. using l a = 1 /-\/ cr /) 
then eqn(§) would hold very accurately. This comparison assumes of course that all 
these gauge theories are similar - which would be true if they were all quite close to 
SU(oo). However we should not assume that this is so at this stage and it would be 
reassuring to have some explicit evidence for the assumption of linear confinement, and 
the universal Luscher correction, for SU(4) and SU(5). Accordingly we show in Fig.|I| 
how the flux loop mass varies with its length in the SU(4) lattice gauge theory, for 
(3 = 10.70. We see some explicit evidence that for L > 10 eqn(Q) is indeed valid, for 
both k = 1 and k = 2 flux strings. For smaller L the higher order corrections to this 
relation are significant, especially so for the k = 2 string whose string tension appears 
to approach the fundamental string tension as its length shortens. We have performed 
a separate calculation (although only for k = 1 strings) at the smaller lattice spacing 
corresponding to (3 = 10.9, and again, as we see in Fig]|, the relation in eqn(fj) holds 
well for the volumes used in Table |T[ 
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With the standard plaquette action that we are using, the leading lattice correction 
to the continuum limit is 0(a 2 ), and so for sufficiently small a we can extrapolate to 
the continuum using 

a k {a) a k (0) 2 

In Fig|| we plot our calculated values of a 2 cjfc = g/a 2 cr/ = cr fc=2 /<7/ against the simulta- 
neously calculated value of aV/. On such a plot the continuum extrapolation in eqn(|3|) 
is a simple straight line. It is clear that all our points are consistent with lying on such 
a line and we show the best fit in the plot. From this we extract our continuum limit 
for the ratio of string tensions: 

a k= 2 _ j 1.354 ±0.070 (3 > 10.55 
a™ a f \ 1.405 ±0.080 (3 > 10.70 ' [ ) ' [ ) 

Both of the fits in eqn(^) are quite good, with confidence levels ~ 40%. (The error 
limits correspond to confidence levels ~ 10%.) The (3 > 10.70 fit has the virtue of 
avoiding the f3 = 10.55 calculation which might be biased by the nearby bulk phase 
transition, but suffers from a shorter lever arm. 

Our SU(5) calculations, which are more limited in scope, are summarised in Table |2] 
and the results for a k =2/o'f are plotted in Fig.|3]. The continuum extrapolation, using 
eqn(^|), is clearly less constrained than for SU(4) and leads to: 

Km = 1.556 ± 0.100 : SU(5). (5) 

In addition we have deliberately tuned the SU(5) couplings so that each of the three 
values of a roughly coincides, when expressed in physical units (of 1/a/ct/), with one 
of the SU(4) couplings. It is plausible that the lattice corrections largely cancel in the 
relative SU(4) and SU(5) ratios, so that the average 

p^/g/W) = i.io ±0,03 : 0.38 > a^r f > 0.24 (6) 

should be close to its continuum value. 



3 Discussion 

According to the MQCD conjecture in eqn(|l]) we should expect 

r o su(2) 
^i-^S-2cos--J 1 su(3) m 

<r f sinf ~ ZC ° S N 1.41... SU(4) Uj 

{ 1.61... SU(5) 
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This is plausible for SU(2), where two fundamental charges can form an unconfmed 
colour singlet, and for SU(3), where two fundamental strings can couple to a single 
fundamental (anti-)string. (We say 'plausible' because one could imagine the N > 4 
string bound states continuing down to unstable, resonant string states for N = 3 and 
N = 2.) Our SU(4) results are consistent with eqnfl7|), within fairly small errors. So are 
our SU(5) results; although the continuum extrapolation is less accurate. The ratio in 
eqn(||) is indeed close to the MQCD prediction of ~ 1.14. All this appears to provide 
significant support for the approximate validity of the MQCD conjecture. 

To gauge how significant this apparent agreement is, it would be useful to compare 
with something else. Since our numerical calculations are on a Euclidean lattice with a 
plaquette action, one obvious comparison is with the strong coupling prediction. This 
is well-known to be 



SU(2) 
SU(3) 
SU(N > 4) 



— = min{fc, N - k} k = { 1 SU(3) (8) 
°7 



which starts out well but is quite inconsistent with our SU(4) and SU(5) calculations. A 
more intriguing comparison, emphasised in M , is with the Hamiltonian strong coupling 
result 

r o su(2) 







SU(2) 


a k k(N - k) k=z 




SU(3) 


a f N-l 


1.33 


SU(4) 


k 1.50 


SU(5) 



This is clearly compatible with our numerical results. Of course one might question the 
relevance of comparing Euclidean lattice calculations with the strong coupling limit of 
the Hamiltonian lattice approach. After all the strong coupling limit is not universal 
and we could find almost anything we wanted in some strong coupling limit or other. 
(For example, if we were to construct a lattice action that uses link matrices in matrix 
representations other than just the fundamental.) On the other hand the origin of the 
formula in eqn(^) has a nice simple physical interpretation that might arise in a more 
general context, and one might argue that colour-electric flux is better represented 
with time continuous. It is therefore particularly interesting to note that our errors 
are already comparable to the difference between the MQCD and Hamiltonian strong 
coupling predictions. A calculation that is better than ours by, say, a factor of O(10) 
would readily discriminate between eqn(|§|) and eqn(^). Such a calculation is quite 
feasible and we intend to carry it out in the near future. 

In any case, we have shown that bound k = 2 strings states do exist for SU(iV > 4) 
and that the binding is strong for SU(4) and SU(5). This will have implications for 
the confinement dynamics. For example [], in the picture where confinement in SU(iV) 



gauge theories is driven by the condensation of centre vortices |Tj| the value of ak=2/&f 



1 this argument arose in a discussion with C. Korthals Altes 
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provides information on the relative density oik = 2 and k — 1 vortices. If one neglects 
all correlations between vortices, then the calculation of how vortices disorder Wilson 
loops involves only a simple Poisson distribution [IS] and one readily finds that the 
number of k = 2 vortices piercing unit area (of the minimal surface spanning a Wilson 
loop) is related to the number of k — 1 vortices by 



Nt =2 _ a k=1 1 



(10) 



N v k=1 a k=2 2 

in the case of SU(4). So using our calculated string tensions we obtain a constraint on 
the relative density of doubly charged centre vortices condensed in the SU(4) vacuum: 
N2/NI ~ 0.25(5). It is interesting to note that at high temperatures one can calculate 
the dual 't Hooft string tension in perturbation theory, and one finds |20 a formula 
that happens to be identical to the Hamiltonian strong coupling result given in eqn(|9]). 

The fact that the k = 2 string tension is much less than twice the fundamental string 
tension, for N = 4 and N = 5, suggests the possibility that the mass spectra of these 
gauge theories might contain prominent states based upon closed loops of such k = 2 
strings, that are additional to any states based on closed loops of fundamental flux. 
It would clearly be interesting to perform lattice calculations of SU(iV) gauge theories 
that are designed to explore this intriguing possibility. 
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SU(4) 


P 


lattice 


MC sweeps 






10.55 


8 4 


2 x 10 5 


0.3726(30) 


1.432(35) 


10.70 


10 4 


10 5 


0.3070(12) 


1.380(19) 


10.70 


12 4 


10 5 


0.3058(12) 


1.377(35) 


10.90 


12 4 


10 5 


0.2421(17) 


1.358(42) 


11.10 


16 4 


9.6 x 10 4 


0.2022(11) 


1.413(42) 



Table 1: The ratio of k = 2 and k = 1 string tensions in SU(4) from lattice calculations 
with various lattice spacings (given in units of the fundamental string tension). 



SU(5) 


P 


lattice 


MC sweeps 






16.755 


8 4 


10 5 


0.3845(20) 


1.54(6) 


16.975 


10 4 


1.6 x 10 5 


0.3027(21) 


1.49(7) 


17.270 


12 4 


8 x 10 4 


0.2480(11) 


1.56(5) 



Table 2: The ratio of k = 2 and k = 1 string tensions in SU(5) for lattice calculations 
with various lattice spacings (given in units of the fundamental string tension). 




Figure 1: The masses of the fundamental, •, and k — 2, o, flux loops versus their 
length, at (3 = 10.7 for SU(4). Shown also is the variation one expects from eqn(Q). 
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m k=1 




L 

Figure 2: The mass of the fundamental flux loop as a function of its length, at (3 
for SU(4). Shown also is the variation one expects from eqn(0). 
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a 



k=2 



a 



f 




a 2 a f 

Figure 3: The ratio of k = 2 and k = 1 string tensions in our SU(4)(«) and SU(5) 
(o) lattice calculations plotted as a function of a 2 <jf. Extrapolations to the continuum 
limit, using a leading 0(a 2 ) correction, are displayed. 
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